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Abstract: Motivated by hints for a light Standard Model-hke Higgs boson and a shift 
■ in experimental attention towards electroweak supersymmetry particle production at the 

CERN LHC, we update in this paper our precision predictions at next-to-leading order of 
perturbative QCD matched to resummation at the next-to-leading logarithmic accuracy for 
direct gaugino pair production in proton-proton collisions with a center-of-mass energy of 
8 TeV. Tables of total cross sections are presented together with the corresponding scale 
and parton density uncertainties for benchmark points adopted recently by the experimen- 
tal collaborations, and figures are presented for up-to-date model lines attached to them. 
Since the experimental analyses are currently obtained with parton showers matched to 
multi-parton matrix elements, we also analyze the precision of this procedure by compar- 
ing invariant-mass and transverse-momentum distributions obtained in this way to those 
obtained with threshold and transverse-momentum resummation. 
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1 Introduction 

For many years, the Standard Model (SM) of particle physics has passed all experimental 
tests. Only the mechanism of electroweak symmetry breaking and of the generation of mass 
has remained as an unsolved question and is now being addressed with the general-purpose 
experiments ATLAS and CMS at the CERN LHC. The discovery of a Standard Model-like 
scalar Higgs boson of mass of about 125 GeV and first measurements of its decay channels 
seem now imminent [1, 2], which would represent an impressive success of this large technical 
endeavor. However, the mass of a fundamental scalar particle is affected by large quantum 
corrections, so that immediately the question of its stabilization with respect to the Planck 
scale arises. 

Weak-scale supersymmetry (SUSY) [3, 4] has long been known to effectively solve not 
only this "hierarchy problem", but also a considerable number of other problems such as the 
Grand Unification of gauge symmetries and the question of a viable dark matter candidate. 
LHC searches for supersymmetric particles have concentrated up to now on the largest, 
strong production channels, but so far no signs of squarks and gluinos have been detected 
[5]. Based on LHC searches with about 1 fb~^ of data, benchmark points have therefore been 
defined for future SUSY searches at the LHC [6], and the strong production cross sections 
and their theoretical uncertainties have been updated to next-to-leading order (NLO) and 
next-to-leading logarithmic (NLL) accuracy for pp collisions with a center-of-mass energy 
of 7 TeV [7]. 

In 2012, the LHC center-of-mass energy has increased to 8 TeV. In addition, the negative 
search results for gluino and degenerate squark production and their cascade decays [5] 
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imply that the experimental attention now shifts towards third-generation squarks and 
the weak production channels, in particular the direct production of gauginos. A first 
investigation of neutralino-chargino pair production and their "golden" decay into three 
charged leptons has already been published by the ATLAS collaboration with the result 
that in simplified models such as the constrained Minimal Supersymmetric Standard Model 
(cMSSM), degenerate lightest chargino (xi) and next-to-lightest neutralino (X2) masses are 
excluded up to 300 GeV, depending on the mass of the lightest SUSY particle (LSP) assumed 
to be the lightest neutralino Xi [8]- Surprisingly, an independent analysis of ATLAS and 
CMS measurements of the W~^W~ cross section seems to fit theoretical predictions better 
when electroweak gauginos of about 100 GeV are included than by the SM cross section 
alone [9]. 

The observations above motivate us to update in this paper our precision (NLO+NLL) 
predictions for direct gaugino pair production. They have been obtained by combining a 
leading order (LO) calculation [10] and its next-to-leading order corrections [11], which had 
already been computed independently previously [12], with resummation of the leading 
and next-to- leading logarithms to all orders in the threshold [11] and small transverse- 
momentum (pt) [13] regimes or simultaneously in both [14]. A similar program has been 
carried out for the other weak channel, i.e. the production of sleptons [15-18]. However, 
their production cross sections are considerably smaller and are not yet accessible at the 
LHC, so that we leave the corresponding theoretical update for future work. The semiweak 
associated production of gauginos with squarks [19, 20] and gluinos [21-23], which may well 
become important should the energy available at the LHC prove insufficient for the pair 
production of the latter, has so far been computed only up to NLO. 

In Section 2 we present, for a representative selection of benchmark points defined in 
Ref. [6], the total gaugino pair production cross sections and the corresponding scale and 
parton density function (PDF) uncertainties at the current center-of-mass energy of 8 TeV 
in tabular form. For comparison and analyses of the data taken previously at 7 TeV, we list 
the corresponding cross sections and uncertainties in Appendix A. In addition, an overview 
of the absolute size of the total cross sections at 8 TeV, the impact of the higher-order 
corrections, and the theoretical uncertainties can be obtained from the figures shown in 
Section 2 illustrating the cross sections as a function of the physical gaugino mass along 
the model lines attached to the selected benchmark points. 

Since the experimental analyses [8] are currently based on LO Monte Carlo simula- 
tions of the gaugino signal with parton showers [24, 25], matched [26] to hard tree- level 
multi-parton corrections [27] and normalized to the NLO total cross section [12], we take 
the opportunity to analyze in this paper the precision of this procedure. We therefore 
compare in Section 3 invariant-mass and transverse-momentum distributions obtained with 
parton showers only [25] and matched [26] to one or two additional hard jets [27] to those 
obtained with threshold [11] and transverse-momentum resummation [13]. Similar compar- 
isons have been presented previously for the hadroproduction of new neutral gauge (Z') 
bosons and charged Higgs {H^) bosons in association with top quarks in Refs. [28] and 
[29, 30], respectively. 

Finally, we summarize our results and present our conclusions in Section 4. 
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2 Total cross sections at y/s = 8 TeV 

As the ATLAS and CMS experiments are now starting to probe the direct production of 
gaugino pairs decaying, e.g., into three charged leptons [8], we present in this Section an 
illustrative selection of total gaugino cross sections at the current center-of-mass energy of 
^/s = 8 TeV in order to facilitate precise comparisons with the experimental results. The 
corresponding cross sections at ^/s = 7 TeV, relevant for data taken in 2010 and 2011, can 
be found in Appendix A. Further results are available from the authors upon request. 

2.1 Benchmark points 

First SUSY search results at the LHC with an integrated luminosity of about 1 fb~^ of 
data excluded six out of the nine SPS benchmark points that had been in use since about a 
decade [31, 32]. These cMSSM points therefore had to be replaced by new lines [6] where, 
motivated by the constraints derived from the measurements of the anomalous magnetic 
moment of the muon {g — 2)^ and the rare 5—7-57 decay, only models with a positive 
off-diagonal Higgs mixing parameter /i > and tan /3 = 10, Aq = GeV or tan f3 = 40, 
Aq = —500 GeV were chosen. Here, tan/3 denotes the ratio of the vacuum expectation 
values (VEVs) of the neutral components of the two Higgs doublets, and Aq is the universal 
soft trilinear coupling of the Higgs fields to the squarks at the Grand Unification scale. 
In discussions among the SUSY working groups of the ATLAS and CMS experiments and 
with the LHC Physics Center at CERN (LPCC) [33], 49 benchmark points lying on these 
lines were retained for detailed cross section studies [34]. 

To keep the number of figures manageable, we have selected thirteen out of these 49 
points which are included in the two model lines 10.1 and 10.3 with tan /3 = 10 and Aq = 
GeV as well as seven additional points located on a line similar to line 40.1 with tan /3 = 40 
and ^0 = —500 GeV. Concerning the tables with detailed cross section analyses presented 
in Section 2.2 and the comparison with the Monte Carlo predictions investigated in Section 
3, we decide to focus in contrast on one specific benchmark scenario for each of the selected 
lines. The latter are listed in bold face in Tables 1, 2 and 3, where we employ the LPCC 
numbering scheme. 

For the first point, with a low value of tan/3 = 10 and a vanishing universal trilinear 
coupling ^0 = GeV, we have chosen an optimistic scenario with low values for the universal 
scalar mass tuq and for the universal gaugino mass mi/2 at the high scale, the point 1 of 
the LPCC numbering scheme. These values of uiq and mi/2 lead to modest neutralino and 
chargino masses at the electroweak scale of about 150-550 GeV (see Table 4 in Section 
2.2). In Table 1 we also list, together with the universal scalar and gaugino masses, the 
physical gluino and average squark masses, which lie in the 1-2 TeV range and are therefore 
compliant with the current experimental bounds from the ATLAS and CMS experiments. 
The associated production of the lightest chargino with the next-to-lightest neutralino 
X2 leads often to the golden trilepton signature. We therefore also present in Table 1 the 
branching ratios of the second lightest neutralino to first and second generation and (s)tau 
(s)leptons. These quantities indeed allow to deduce the trilepton production rate from the 
total cross sections given in the next subsections, since for all the scenarios lying on the 
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Point 


(mi/2,mo) [GeVJ 


rrig |CjeVJ 


{rriq) [GeVJ 




1)1} ro/1 
11 ITT) [7oJ 


1 


(400, 100) 


935 


840 


16 / 


/ 22 


1 


(45L), ii/.oj 


1040 


940 


20 / 


' 20 


3 


(500, 125) 


1145 


1030 


24/ 


/ 19 


4 


(550, 137.5) 


1255 


1125 


26 / 


/ 18 


5 


(600, 150) 


1355 


1220 


28 / 


/ 18 


6 


(650, 162.5) 


1460 


1310 


28 / 


/ 17 


7 


(700, 175) 


1565 


1405 


29 / 


/ 17 



Table 1. Selection of benchmark points on the cMSSM model line 10.1 of Ref. [6] with tan /3 = 10, 
^0=0 GeV, and toq = 0.25 x mi/2- The points are spaced in steps of Ami/2 = 50 GeV, and the 
gluino and average squark masses are rounded to 5 GeV accuracy. We also present the branching 
ratios of the next-to-lightest neutralino into £ ^ fi and r (s)leptons. 



Point 


(mi/2,mo) [GeV] 


m-g [GeV] 


{mg) [GeV] 


BR(x° ^ x?/i) [%] 


15 


(375, 250) 


885 


825 


89 


16 


(450, 300) 


1050 


975 


92 


17 


(525, 350) 


1210 


1125 


92 


18 


(600, 400) 


1370 


1275 


92 


19 


(675, 450) 


1525 


1420 


92 


20 


(750, 500) 


1680 


1565 


92 



Table 2. Selection of benchmark points on the cMSSM model line 10.3 of Ref. [6] with tan f3 = 10, 
Aq ~ GeV, and mi/2 — 1.5 x ttiq. The points are spaced in steps of Amp = 50 GeV, and the 
gluino and average squark masses are rounded to 5 GeV accuracy. We also present the branching 
ratio of the dominant decay mode of the next-to-lightest neutralino into a Higgs boson. 

line 10.1 the lightest chargino and the sleptons always decay (in a first approximation) 
to a single observable lepton and missing transverse energy. 

Our second benchmark scenario consists in the point 18 lying on the model line 10.3. 
As for the point 1, the value of the ratio of the Higgs VEVs has been set to tan /3 = 10 and 
the universal trilinear scalar coupling is vanishing at high energies. In contrast to the first 
scenario, we rather adopt here higher values for the universal scalar and gaugino masses itiq 
and mi/2, which yields slightly heavier neutralino and chargino states with masses of about 
250-770 GeV (see Table 5 in Section 2.2). As for the first scenario, this combination of 
values for mo and mi/2 shifts the masses of the squarks and gluino to the 1-2 TeV range as 
shown in Table 2, rendering the scenarios of the line 10.3 not (yet) excluded by LHC data. 
As shown in the last column of Table 2, the next-to- lightest neutralino decays most of the 
time to a Higgs boson, which is a clear sign of its higgsino nature. Even if its interactions 
with first and second generation squarks and quarks are consequently suppressed, we remind 
that neutralinos (and charginos) couple to the electroweak W and Z bosons through their 
higgsino components. Therefore, one expects reasonable production rates for final states 
containing one or two neutralinos X2 s-channel weak boson exchanges, in contrast to the 



- 4 - 



Point 


(mi/2,mo) [GeV] 


m-g [GeV] 


{rriq) 1 mi^ / m 


I. [GeV] 




27 


(400, 300) 


940 


890 / 615 / 


745 


98/2 


28 


(450, 325) 


1050 


985 / 700 / 


835 


97/3 


29 


(500, 350) 


1155 


1185 / 780 / 


925 


96/4 


30 


(550, 375) 


1260 


1180 / 860 / 


1010 


95/4 


31 


(600, 400) 


1365 


1275 / 940 / 


1100 


95/5 


32 


(650, 425) 


1470 


1370 / 1020 / 


1185 


94/5 


33 


(700, 450) 


1575 


1465 / 1095 / 


1275 


94/6 



Table 3. Selection of benchmark points on a cMSSM model line similar to line 40.1 of Ref. [6] 
with tan/3 = 40, Aq ~ —500 GeV, and mo = 0.5 x + 100 GeV. The points are spaced in steps 
of Atoi/2 = 50 GeV. The mass of the gluino, the average mass of the heavier squarks. and the 
masses of the (considerably lighter) top and bottom squarks are rounded to 5 GeV accuracy. We 
also present the branching ratios for the decays of the next-to-lightest neutralino into r (s)leptons 
and Higgs bosons. 



benchmark point 1 where the second lightest neutralino is mainly a wino and thus rather 
produced through t- and u-channel squark exchanges. Since the lightest chargino decays 
with a 98% branching ratio to a VF-boson, one can deduce the production rate of golden 
signatures from the values of the X2xf total cross section presented in Section 2.2 and 
Section 2.3 after accounting for the branching ratio of a Higgs boson decaying to a pair of 
tau leptons and the rate of leptonic VF-boson decays. 

Our third selected scenario lies on a model line similar to the line 40.1 of Ref. [6] with a 
large value of tan /? = 40 and a large negative value of the universal scalar trilinear coupling 
Aq = —500 GeV. We retain the high values of the universal scalar and gaugino masses 
that have driven the choice of our second scenario and therefore adopt the point 31 of the 
LPCC numbering scheme as our third benchmark point. Comparing with the point 18, the 
neutralino and chargino masses are not drastically affected by the different choices of tan /3 
and ^0 ci'iid hence lie in the 250-825 GeV range (see Table 6 in Section 2.2). Similarly, the 
gluino and first and second generation squark masses are not drastically affected either as 
presented in Table 3. In contrast, the large values of tan/3 and Aq induce an important 
mixing among the third-generation squark interaction eigenstates (and to a smaller extent 
among the left- and right-handed stau states). As a consequence, the masses of the lightest 
stop and sbottom are much lower than the average squark mass as shown in Table 3, 
and these states are almost maximal admixtures of the left-handed and right-handed third 
generation squark eigenstates. We also present the two main decay modes of the second 
neutralino, decaying mainly to an associated pair of stau (always decaying to a tau and a 
lightest neutralino) and tau lepton but also, at a much smaller rate, to a Higgs boson and 
missing energy carried by the lightest neutralino. As for the two other selected scenarios, 
the cross section corresponding to the trilepton mode can be deduced from the values of 
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the total cross sections given in Section 2.2 and Section 2.3 and from the branching ratios 
of the hghtest Higgs boson to taus as weh as the one associated to tau leptonic decays. 

We emphasize that while the increasing experimental limits on the squark and gluino 
masses might exclude the corresponding cMSSM scenarios, the (lower) gaugino masses 
may still remain allowed and their total cross sections therefore approximately valid in the 
context of more general SUSY-breaking scenarios. 

2.2 Total cross sections for the benchmark points 1, 18, and 31 

The production of pairs of charginos and neutralinos has been initially studied at leading 
order of perturbative QCD in the early 1980s [35, 36], while more recently polarization 
[10, 37] and flavor-violating [38-40] effects have been included and investigated. Next- 
to-leading order corrections have been supplemented to the leading order approximation 
[11, 12] and have been found to be important, in particular due to the presence of large 
logarithmic contributions arising from soft and collinear parton emission by the initial state 
particles. Since they spoil the convergence of the perturbative series, these logarithms have 
to be resummed to all orders in the strong coupling constant to allow for reliable theoretical 
predictions in the entire phase space, including the regions related to soft and collinear QCD 
radiation. Transverse-momentum and threshold resummation have then been achieved at 
the next-to-leading logarithmic accuracy independently [11, 13] and simultaneously [14]. 

In this Section, we present in Tables 4, 5, and 6 total production cross sections for 
various combinations of neutralino and chargino pairs in the context of the benchmark 
points 1, 18 and 31 introduced in Section 2.1 at the leading and next-to- leading order of 
perturbative QCD as well as after matching the NLO results to threshold resummation. 
The masses of the produced gauginos are also given as references, while the values of the 
squark masses entering the partonic cross section can be found in Tables 1, 2 and 3. 

The LO results, given in the fourth column of the Tables, are computed as in Ref. 
[10] after convolving (unpolarized) partonic cross sections with the LO set of the MSTW 
2008 parton densities [41], accounting for five light flavors of massless quarks and as agreed 
among the SUSY working groups of ATLAS, CMS and the LPCC. We use a top quark 
mass of 173.1 GeV [42] and the values of mz = 91.1876 GeV and mw = 80.403 GeV [43] 
for the masses of the electroweak gauge bosons, and we consider the CKM matrix as the 
identity matrix in flavor space. The supersymmetric spectra have been generated with 
the SuSpect 2.41 program [44], allowing to obtain low-energy supersymmetric masses and 
parameters from universal parameters defined at the Grand Unification scale and evolved 
down through renormalization group running at the two-loop level. The central values of 
the cross sections in the Tables have been computed after fixing the factorization scale to 
the average final state particle mass, while the uncertainties are estimated after multiplying 
the central value of the factorization scale by a factor lying in the 0.5-2 range. 

For our NLO predictions, we follow the detailed computations preformed in our previous 
work [11], where both the QCD and the SUSY-QCD contributions are included in the 
virtual pieces of the cross section, and give our results in the fifth column of the Tables. 
We emphasize that we also consider, as in our earlier calculations, internal squark mixing, 
i.e., the squarks are kept non-degenerate in the loops. The partonic cross sections are 
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Process 



mi [GeV] 



ms [GeVj 



LO [fb] 



NLO [fb] 



NLO+NLL [fb] 



PP 
PP - 
PP 
PP 
PP - 
PP 
PP - 
PP - 
PP - 
PP - 
PP - 
PP - 
PP - 
PP - 
PP - 
PP - 



• XiXi 

>X'2XI 

■ X2X1 

^ A3A4 

■ X3X1 

■ X3X2 

■ X4X2 

■ xtxi 

■ xtX2 

■ xtxs 

xtxT 

■ xtXs 

■ xtxl 

xtx2 



161.7 
161.7 
303.8 
303.8 
303.8 
526.5 
526.5 
526.5 
542.4 
303.5 
303.5 
303.5 
303.5 
542.2 
542.2 
542.2 



161.7 
303.5 
303.8 
526.5 
303.5 
542.4 
303.5 
542.2 
542.2 
161.7 
303.8 
526.5 
303.5 
526.5 
542.4 
542.2 



0-8115:3% 
^•-^"-5.5% 

n 91+9.4% 

^•^-^-8.1% 

14 46+6-7% 
n sQ+ii-0% 

U-0'3_9.3% 
r. 1 9+9.4% 

Q 42+11.2% 

^•^^-9.5% 
Q on+11.3% 
^•■^^-9.6% 
n qo+6.0% 
0-='«-5.4% 

34+9-2% 

9^ fi4+6-6% 
^0-04_5.9% 

-1 97+11.1% 
-^•^'-9.4% 
-1 9-1+11.2% 
^•^-^-9.5% 

0-80-9.3% 



1.06 

0.20 
1.07 
0.25 
17.25 
0.97 
0.15 
0.50 
0.47 
0.46 
40.90 
0.40 
30.37 
1.46 
1.37 
1.00 



+3.5% +2.8% 
3.0% -2.0% 
+2.5% +2.9% 
2.4% -2.4% 

+3.5%+3.1% 
3.5% -2.2% 
+2.6% +3.2% 
2.9% -2.3% 
+1.6%+3.0% 
1.7% -2.6% 
+2.8% +3.9% 
3.3% -2.4% 
+2.6% +3.8% 
2.9% -2.9% 
+2.8% +4.9% 
3.3% -3.6% 
+2.7% +4.9% 
3.2% -3.6% 
+2.5% +2.8% 
2.4% -2.1% 
+1.6%+2.9% 
1.7% -2.2% 
+2.6% +3.7% 
2.9% -2.4% 
+1.7%+2.7% 
1.9% -2.0% 
+2.9% +4.4% 
3.3% -2.7% 
+2.7% +4.4% 
3.2% -2.8% 
+2.6% +4.0% 
3.1% -2.4% 



1.03 
0.20 
1.05 
0.25 



+0.5% +2.9% 
-0.6% -2.0% 
+0.0% +2.9% 
-0.3% -2.5% 

+0.0% +3.5% 
-0.4% -1.9% 
+0.1% +3.2% 



17.05 
0.96 
0.15 
0.49 
0.46 
0.46 



-0.5% -2.3% 

+0.2%+3.1% 
-0.7% -2.6% 
+0.4% +3.8% 
-0.9% -2.5% 
+0.1% +3.8% 
-0.6% -3.0% 
+0.4% +4.9% 
-0.9% -3.5% 
+0.5% +4.9% 
-l.l%-3.7% 
+0.2% +2.9% 
-0.5% -2.1% 



40.51 



+0.0% +2.9% 
-0.3% -2.2% 



0.40 



+0.0% +3.6% 
-0.3% -2.5% 



30.04 



+0.0% +2.7% 
-0.5% -2.1% 



1.45 
1.36 
0.99 



+0.3% +4.3% 
-0.7% -2.9% 

+0.4% +4.6% 
-0.8% -2.6% 
+0.4%+4.1% 
-0.9% -2.4% 



Table 4. Total cross sections related to the production of various gaugino pairs of masses mi 
and m2, presented together with the associated scale and PDF uncertainties for the LHC running 
at a center-of-mass energy of y/s = 8 TeV in the context of the benchmark point 1 of the LPCC 
numbering scheme. The cross sections are given at the leading order and next-to-leading order of 
perturbative QCD and matched to threshold resummation. The PDF uncertainties are not shown 
for the LO results. Any cross section smaller than 0.1 fb is omitted. 



Process 


mi [GeV] 


7712 [GeV] 


LO [fb] 


NLO [fb] 


NLO+NLL [fb] 


PP XiXi 
PP X2X1 
PP xtX2 

PP xtxi 

PP xtXs 

PP xtxl 

PP xtX2 


249.6 
471.9 
471.8 
471.8 
766.3 
766.3 
766.3 


249.6 
471.8 
471.9 
471.8 
754.0 
766.6 
766.3 


n 1 q+8.6% 

4 70+9.8% 

^•'"^-8.4% 
q 1 q+9.8% 
>3-J--3_8.4% 
Q ig+14.2% 
^•■^"-11.6% 
n 1 t:+14.3% 

n 1 1 +13.6% 
^•^^-11.2% 


n iR+3.5%+3.3% 
^•-^"-3.4% -2.3% 
1 OQ+1-8%+4.1% 
^•°°-2.4%-3.1% 
c 9o+1.8%+3.9% 
d.ZO__2.4%-2.5% 
q c7+1.9%+3.5% 
■^•^'-2.5%-2.2% 
n i7+3.5%+6.1% 
"•■^ -4.2% -3.8% 

n iR+3.4%+6.1% 
U.iD_4 2%_3.9% 

n 19+3.1% +6.0% 
"•^^-3.9% -3.5% 


n ifi+0.2%+3.5% 

^•-^"-0.3% -2.4% 

1 of;+0.6%+4.1% 
^•°"-1.2%-3.1% 
c 99+0.3% +4.0% 
'J-^^_0.6%-2.5% 
q c;9+0.4%+3.7% 
•^•^■^-0.7% -2.3% 
n i7+1.0%+6.1% 
"•-^ -1.8% -3.8% 

n iR+i.1%+6.1% 

^•-^"-1.8% -4.0% 
n i9+1.0%+6.0% 
"•-^^-1.8% -3.6% 



Table 5. Same as Table 4 for the benchmark point 18 of the LPCC numbering scheme. 



convolved this time with the MSTW 2008 NLO parton density sets [41], and the central 
values of the total cross sections are given in the Tables together with the theoretical 



Process 


nil [GeV] 


1712 [GeV] 


LO [fb] 


NLO [fb] 


NLO+NLL [fb] 
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n. 1 c-+0.2%+3.4% 
^•-^^-0.3% -2.4% 
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^•O-L_0.6%-2.6% 
q 94+0.5% +3.8% 
■^•^^-0.7% -2.3% 



Table 6. Same as Table 4 for the benchmark point 31 of the LPCC numbering scheme. 



uncertainties related to scale variation and the choice of the PDF set. Scale uncertainties 
are derived following the standard approach of simultaneously varying the factorization 
scale and renormalization scale fiji by multiplying and dividing the average final-state 
particle mass by a factor of two. The PDF errors are obtained by evaluating the envelope 
of the cross section when employing the 68% confidence level range of the MSTW 2008 
parton densities. The (asymmetric) errors AcJup and Acidown are computed as defined by 
the MSTW collaboration, i.e., according to 

n 2 

(AcJup)^ = ^ |max[cr^ - ctq, - ctq, O] | , 

2 

(AtTdown)^ = |max[o-o - ao-cr^, O] | , 

k=l 

where (Tq consists in the value of the cross section computed when using the central set of 
parton densities, while and are those obtained from ztlo" variations along the /c*^ 
eigenvector of the covariance matrix associated to the PDF fit. 

In the last column of the Tables, we present the total cross section after having matched 
the next-to-leading order results with threshold resummation at the next-to-leading loga- 
rithmic accuracy. After performing a Mellin transformation, the cross section can be written 
in the Mellin A^-space as a simple product between the parton densities /^/p and f^/p and 
the partonic cross section aab, 

1) = ^fa/p{N,fi^)fb/piN,fi^) aab{N,M\i?) , (2.2) 

ab 

where we have set /i_p = fiR = /i for brevity and where M denotes the invariant mass of the 
produced gaugino pair. Within the threshold resummation formalism, the partonic cross 
section can be refactorized into a closed exponential form as [45] 

aab{N,M^,fi^)=nab{M^,fi^) exp [OabiN, , fi^)\ • (2.3) 

The function Ti can be perturbatively computed, describes the hard part of the scattering 
process, and is independent of the Mellin variable N. The function Q collects soft and 
collinear parton emission and absorbs the large logarithmic contributions arising at fixed 
order. We then further improve the resummation formula of Eq. (2.3) above by including 
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in the ?^-function the dominant 1/A^-terms stemming from universal cohinear radiation of 
the initial state partons [46-49] . The exact expressions of the coefficients of the expansion 
of the two functions and Q at the next-to-leading logarithmic accuracy can be found in 
Ref. [11]. 

The resummed predictions are known to be valid near production threshold, where the 
logarithmic terms of the cross section dominate, in contrast to the fixed order computations, 
which are assumed to be valid far from this threshold. Therefore, in order to get reliable 
predictions in all kinematic regions, both results are consistently matched by summing the 
esumi 

(exp.) 



resummed a^^^'^ and the fixed order cr^^j,'" '' cross sections and subtracting their overlap 



(res.) . (f.o.) (exp.) /„ ^\ 

^ab = <fe ' + Kb ' - <b ' ■ (2.4) 

We again refer to Ref. [11] for the exact expression of o'^j,^^ '* which is, for a NLO+NLL 
matching, the 0{as) series expansion of the resummed cross section of Eq. (2.3). Since the 
two quantities a'^^^'^ and c^^^''''' are computed in Mellin space, an inverse transform must 
be performed in order to obtain the total cross section in terms of physical quantities 

^2_da_^_^ _ 1 /■ ^j,r^~Nf fAT , ..2\f /Tvr , 1 ,.2^ 

d^^' ---■'c. (2.5) 



M^^(r) = 1^Y.I dNT-^'fa/piN + 1, f^^)h/piN + 1, fi 

ab ''^f^ 



where r = M"^ /s. To avoid singularities in the integrand related to the Landau pole of 
the strong coupling constant and to the Mellin moments of the parton densities at small 
momentum fraction, we choose an integration contour Cm according to the principal value 
procedure [50] and the minimal prescription [51]. 

The resummed results are given in Tables 4, 5 and 6 together with the uncertainties 
derived from scale variations and deviations from the central fit of the MSTW 2008 parton 
densities. The latter are computed as for the pure NLO results and described above. 

One observes, by investigating the numerical results, an improvement of the stability 
of the perturbative series, as the scale dependence of the total cross section is tamed by 
the resummation of the soft and collinear radiation contributions. At LO, the factorization 
scale fjip appearing within the evolution of the parton densities already induces large loga- 
rithmic terms yielding an uncertainty of about 10% on the cross section. This effect is then 
attenuated at NLO, but one gets an additional explicit dependence in the renormalization 
scale ///J through the strong coupling constant and the loop contributions. After matching 
the NLO results with threshold resummation at NLL, the scale variations are eventually 
considerably reduced thanks to the inclusion of dominant higher order contributions within 
the Sudakov form factor Q. In contrast, the uncertainties related to the choice of the 
employed parton densities coincide when comparing the NLO and the NLO+NLL results, 
which is not surprising since the same PDF sets enter both computations. 

2.3 Total cross sections for model lines 10.1, 10.3, and 40.1 

In order to give a more complete overview of the absolute size of the total gaugino cross 
sections, we present in this Section the total cross sections for the golden trilepton channel. 
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300 350 400 450 500 550 

m~ [GeV] 

Figure 1. The total cross section at LO (dotted), NLO (dashed), and NLO+NLL (fuU) with its 
scale (green) and PDF (yellow) uncertainty for the production of a X2X1 P^-ir ^ function of their 
almost degenerate mass rrij^ at the LHC with \/s = 8 TeV and benchmark points 1-7. 



i.e. the production of a paii", as a function of their almost equal mass m-j^ at the LHC 

with a center-of-mass energy of 8 TeV 

The results for model line 10.1 are shown in Figure 1 at LO (dotted), NLO (dashed), 
and NLL+NLO (full). As the soft SUSY- breaking mass mi/2 increases from 400 to 700 
GeV, the gaugino masses rise from 300 to 550 GeV. Consequently, the total cross section 
falls by over an order of magnitude from about 40 fb to about 2 fb. The impact of the NLO 
corrections to the LO prediction is clearly visible, and the cross section is then stabilized 
at NLL+NLO, where it remains almost unchanged with respect to the NLO result, but is 
considerably less scale-dependent. The overall theoretical uncertainty is thus dominated by 
the PDF uncertainty, estimated as described in the previous Section. 

Similar results for model lines 10.3 and 40.1 can be seen in Figure 2 and Figure 3, 
respectively. Despite a very similar physical gaugino mass range, the total cross sections 
are larger (falling from about 60 to 3 fb) for model line 10.3 than for line 10.1, which can 
be attributed to the higgsino nature of the sparticles, so that they are now also produced 
through s-channel weak gauge bosons. On the other hand, the total cross sections for model 
line 40.1 are very similar (falling from about 35 to 2 fb) to those of line 10.1. The size of the 
corrections and the uncertainties show no substantial differences among the three scenarios. 
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P P ^ X % the LHC (8 TeV), scenarios 15-20 




300 350 400 450 500 

m~ [GeV] 

Figure 2. Same as Figure 1 at benchmark points 15-20. 

3 Comparisons of invariant-mass and transverse-momentum distribu- 
tions 

Typical final-state signatures associated to the production of the electroweak superpartners 
of the Standard Model gauge and Higgs bosons contain in general abundant initial state 
QCD radiation. The latter has important effects on the shapes of the kinematical distri- 
butions, in particular due to the large logarithmic contributions arising in the phase space 
regions where these additional partons are neither widely separated nor hard. In this case, 
reliable theoretical predictions require a consistent reorganization of the logarithmic terms 
so that they are effectively resummed to all orders in the strong coupling constant and em- 
bedded in the so-called Sudakov form factor. In the resummation computations presented 
in Section 2, this factor has been derived at the next-to- leading logarithmic accuracy from 
the knowledge of the perturbative expansion of the hard-scattering function at the next- 
to-leading order. This leads to an accurate description of the QCD environment in phase 
space regions where radiation is soft and/or collinear. In order to obtain a good description 
over the whole kinematical range, including also the hard radiation regime, the resummed 
results have to be eventually matched to a fixed order calculation describing precisely the 
hard emission. 

In experimental analyses or phenomenological investigations relying on Monte Carlo 
simulations, the production of additional jets is traditionally simulated using parton show- 
ering programs, such as Herwig [24, 52] or Pythia [25, 53], which describe QCD emis- 
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Figure 3. Same as Figure 1 at benchmark points 27-33. 



sions as successive branchings of a mother parton into two daughter partons. These tools 
are based on Markov chain techniques built upon the Sudakov form factor dictating the 
probability laws for a specific parton to radiate or not. In contrast to the resummation 
techniques adopted for the predictions shown in Section 2, where the Sudakov form factor 
is computed at the next-to-leading logarithmic accuracy, the parton showering tools use 
only a leading-logarithmic precision, sometimes including partial (but not complete) next- 
to-leading logarithmic contributions. The reason is that only such approximations allow for 
a description in terms of a Monte Carlo algorithm. Moreover, for supersymmetric processes 
the hard scattering can only be evaluated at tree level, since no SUSY process has so far 
been implemented into NLO Monte Carlo tools such as MC@NLO [54, 55] or POWHEG 
[56]. This description is formally only correct in the soft and collinear radiation regions of 
the phase space, and it fails when considering the production of hard and widely separated 
additional partons due to missing subleading terms. In this case, matrix elements describing 
the same final state together with an additional parton are required. 

The matrix-element approach is indeed accurate for describing hard and widely sep- 
arated emissions, but is known to break down in the soft and collinear radiation limits. 
Contrary, parton showering gives a proper description in the soft and collinear kinematical 
regions, but underestimates hard emissions. Therefore, several matching algorithms, such 
as the Catani-Krauss-Kuhn- Webber (CKKW) scheme based on event reweighting [57, 58] 
or the Maganano (MLM) scheme based on event rejection [59], have been developed to 



combine the two methods in a consistent way. This hence avoids the possible double count- 
ing arising from radiations which can be taken into account both at the level of the matrix 
element and at the level of the Sudakov form factor. These different procedures have been 
extensively applied, compared and confronted to the data in the context of various Standard 
Model processes [60-63] and also used for recent Beyond the Standard Model explorations 
[64-66]. 

In this Section, we present the first comparison between resummation (matched to fixed 
order), next-to-leading order and parton showering (matched to matrix elements) predic- 
tions in the context of chargino and neutralino production. We focus on the distributions 
of the transverse momentum (px) and invariant mass (M) of the superparticle pairs pro- 
duced at the LHC collider, running at a center-of-mass energy of 8 TeV. Concerning the 
most precise predictions, we employ the transverse- momentum and threshold resummation 
formalisms. Threshold resummation has been briefly described in Section 2.2 and more 
details can be found in Ref. [11]. As for the threshold case (see Eq. (2.2)), resummation in 
the transverse-momentum regime is also performed in Mellin space. The partonic cross is 
computed in the impact parameter 6-space 

aabiPT, M\fl^) = / db-Jo{bpT) CTabib, N, M\ S, I?) , (3.1) 

where the function Jq is the 0*'^-order Bessel function, so that the integrand can be refac- 
torized under a closed exponential form, as in Eq. (2.3) for threshold resummation, 

aab{h,N,M^,^?)=nab{N,M\^Ji^) exp kf,(6, iV, M^, /x^)! . (3.2) 



As usual, the Sudakov form factor Q contains the soft and collinear radiation contributions 
and 7i is the hard function, independent of the impact parameter b. More details and the 
exact expressions of these functions, evaluated at the next-to-leading logarithmic accuracy, 
can be found in Ref. [13]. Matching to the fixed order computations and the inverse Mellin 
transform to get back to the physical space are eventually performed, as sketched in Eqs. 
(2.4) and (2.5). 

Concerning the merging of the parton showering with the hard-scattering matrix ele- 
ments, we follow the fc^-MLM matching scheme [64] as implemented in the MadGraph- 
MadEvent event generator version 5 [27], interfaced to Pythia 6 [25]. The UFO [67] 
model files required for MadGraph are generated starting from the MSSM implementation 
in the FeynRules package [68-70] after loading the supersymmetric spectrum associated 
to the benchmark points 1, 18 and 31 described in Section 2 and obtained as in Section 2.2 
with the SuSpect 2.41 program [44]. In the fc^-MLM matching scheme, matrix-element 
based events are allowed to contain zero, one or two additional partons in addition to the 
pair of supersymmetric particles. These (parton-level) jets are generated with a minimum 
jet measure of 50 GeV between two final state partons i and j, where kx is defined by 

= Ta.m.{p^i,p\j)Rij . (3.3) 

In the equation above, the quantities pxi and pTj are the transverse momenta of the two 
partons under consideration and Rij denotes their angular distance in the (ry, (f)) plane. In 
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the case where one of the two partons is an initial state parton, we define the jet measure 
as the transverse momentum of the final state parton 

kr = PTi (3.4) 

and require its value to be larger than 20 GeV. The events are then passed to Pythia for 
showering, and jets are reconstructed using FastJet [71] in the fey-jet algorithm with a 
cut-off scale Q^'^t'^'^ set to 70 GeV. The jets are said to be matched to one of the original 
partons if the /cr-measure between the jet and the parton is smaller than Q'^^^'^^. The events 
are selected only if each jet is matched to one parton, with the exception of the two-jet 
sample. In this case, extra jets are still allowed, which maintains the full inclusiveness of 
the matched sample. 

In our simulation setup, we neglect all quark masses but the top mass and employ the 
leading order set of the MSTW 2008 parton density fits [41]. In contrast, the NLO and 
NLO+NLL results are based on the next-to-leading order set of the MSTW 2008 fits. We 
identify the factorization and renormalization scales as the average mass of the produced 
particles in the case of the three hard processes, i.e., without any extra parton as well as 
with one and with two additional (hard) radiations, while the renormalization scale used 
by the showering algorithm consists in the jet measure at each branching. The value of the 
total cross section associated to a matched event sample is usually rather close to the one 
of the original process, where no extra radiation is allowed, i.e., it is close to the value of 
the tree-level cross section. Therefore, we reweight the produced events uniformly, so that 
the total rate is now equal to the resummed cross sections given in the Tables of Section 
2. Let us note that we have also performed the on-shell subtraction of any intermediate 
resonance which could appear in the samples containing one or two extra partons using the 
narrow-width approximation, as those subprocesses are rather identified with other genuine 
processes such as the associated production of a gaugino and a squark or a gluino. Finally, 
the event samples have been analyzed with the program package MadAnalysis 5 [72]. 

3.1 Distributions in the invariant mass of the gaugino pair 

In Figure 4, we depict the invariant mass spectrum da/dM for the associated production 
of a lightest chargino (xf) ^^'^ second lightest neutralino (xg) pair in the context of the 
benchmark scenario 1 presented in Section 2.1. We compare the results obtained using 
the threshold resummation formalism (thick red full curve), found to be almost equivalent 
to those corresponding to a pure next-to-leading order calculation, to those generated by 
LO Monte Carlo simulations with MadGraph and interfaced to the parton showering 
algorithm provided by Pythia. For the latter, we have normalized the total cross section 
to its resummed value of 40.51 fb (see Table 4). The spectrum starts at M ~ 600 GeV, since 
the two gauginos are almost degenerate with a mass of about 300 GeV each. The shape of 
the spectrum is then dictated by the nature of the produced superparticles, which are here 
mainly gaugino-like. This suppresses the s-channel VT-boson exchange contributions with 
respect to P-wave production through t- and u-channel squark exchange, the mass of the 
squarks being here of the order of 800 GeV. The resummed contributions increase the total 
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Figure 4. Distributions in the invariant mass M of a X2xt P^i^ with mass 304 GeV each (bench- 
mark point 1) at the LHC with y/s = 8 TeV. We compare the NLO matched to the NLL (red, 
thick full) distribution to the results obtained after matching matrix elements containing no (green, 
dotted), one (blue, dashed), and up to two (red, dot-dashed) additional jets to parton showering. 

rate in comparison to the pure next-to-leading order results only slightly, since threshold 
effects only become important when the invariant mass approaches the total center-of-mass 
energy. 

The set of three curves obtained using the MadGraph and Pythia programs agree 
quite well with each other. They are related to three different ways which we have adopted 
to perform the matching procedure. The results represented by the green dotted curve have 
been obtained from tree-level matrix elements only, related hence to the process 

PP X2xt (3-5) 

without any extra radiation. The parton-level events have then been passed to Pythia 
for parton showering, and no matching procedure has been applied. After accounting for a 
resummed X-factor normalizing the integrated distribution to 40.51 fb (see Table 4), one 
observes a good agreement with the resummed predictions. A similar behavior has already 
been observed in the context of Z' production at the LHC, where in Ref. [28] resummed 
computations were confronted to results obtained with tree-level matrix-elements passed to 
a parton showering algorithm. 

Applying the /ct-MLM matching procedure described in the beginning of Section 3, we 
compute matrix elements containing up to one additional QCD emission, i.e., considering 
the subprocesses 

PP X2xt and pp ^2 xfj^ (3-6) 
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Figure 5. Same as Figure 4 for the benchmark point 18. The gaugino masses are here 472 GeV. 

where j denotes a quark, an antiquark or a gluon, and match them consistently to par- 
ton showering. We obtain the results shown by the blue dashed curve of Figure 4 after 
normalizing again the total cross section to its resummed value of 40.51 fb (see Table 4). 

In a similar fashion, one can allow for matrix elements with up to two extra partons, 
i.e., we account for the processes 



This leads to the dot-dashed red curve presented in Figure 4. Again, since the matched 
total cross section is close to the unmatched one, the results have been normalized to 40.51 
fb, including hence a resummed X-factor. 

In the case of the benchmark scenarios 18 and 31, the gaugino masses are heavier, with 
rriy. = 472 GeV and 479 GeV, respectively. This reduces the total production cross section 
of an associated pair of the lightest chargino and the next-to-lightest neutralino by about 
an order of magnitude, which then reads 5.22 fb for the scenario 18 (see Table 5) and 4.81 
fb for the scenario 31 (see Table 6). Employing these values for the normalization of the 
results generated by MadGraph and Pythia, we confront them again to our NLL+NLO 
calculation in Figures 5 and 6. One observes a similar behavior as for the scenario 1 shown 
in Figure 4. 

3.2 Distributions in the transverse momentum of the gaugino pair 

For scenario 1 at the LHC with a center-of-mass energy of 8 TeV, we show in Figure 
7 the transverse- momentum distribution da/dpT of a lightest chargino (xt) and next- 



PP X2 Xi , PP^ X2 Xi i, and pp X2 Xi J j ■ 



(3.7) 
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Figure 6. Same as Figure 4 for the benchmark point 31. The gaugino masses are here 479 GeV. 

to-lightest neutralino (xg) pair- As expected, the fixed order predictions at 0{as) (full 
blue curve) diverge as the transverse momentum pT tends to zero due to the unbalanced 
large logarithmic terms related to soft parton radiation. After their consistent matching 
to a resummation computation in the transverse-momentum regime at the next-to-leading 
logarithmic accuracy (thick red full curve), the results exhibit a finite behavior with a 
maximum around a p^^-value of about 10 GeV. The effects of the resummation of the 
large logarithms extend to intermediate values of the transverse momentum, where the 
resummed predictions are considerably larger than those obtained from the pure fixed order 
computation as it has already been observed for slepton-pair production [16, 18] or other 
gaugino-pair production channels [13, 14]. 

We now turn to the confrontation of the resummed results to the three sets of predic- 
tions obtained after employing the MadGraph and Pythia generators, following the same 
approach as the one described in Section 3.1. Again the total production rate associated 
to each of these three curves has been set to the value of the resummed total cross section 
of 40.51 fb (see Table 4). One observes very good agreement between the most precise 
distribution including a matching of the resummed results to the fixed order ones (thick 
red full curve) and the distributions obtained with MadGraph and Pythia after having 
matched the matrix elements containing one or two additional jets to parton showering 
(blue dashed and red dot-dashed curves in Figure 7) . Contrary, the simple application of a 
parton showering algorithm to the tree-level matrix element (green dotted curve on Figure 
7) leads to a slightly too soft spectrum as expected, since parton showering methods can 
not address properly the intermediate and large pT regions. 
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Figure 7. Distributions in the transverse momentum pT of a X2X1 P^'ii' with mass 304 GeV each 
(benchmark point 1) at the LHC with ^/s = 8 TeV. We compare fixed order at 0{as) (blue, full) and 
NLL (red, thick full) distributions to the results obtained after matching matrix-elements containing 
no (green, dotted), one (blue, dashed), and up to two (red, dot-dashed) additional jets to parton 
showering in the small (top) and large (bottom) px regions. 
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In contrast to the invariant-mass distributions, where the three Monte Carlo predictions 
(without matching, with a matching procedure of matrix elements accounting for at most 
one additional parton, and with a matching procedure of matrix elements containing at 
most two extra partons) have been found to be in a rather good agreement (see Section 
3.1), one can clearly here observe the effects of including additional partons when comparing 
the three distributions generated by the Monte Carlo simulators in Figure 7. While the 
position of the peak does not depend much on the presence of additional partons at the 
matrix-element level, the global hardness of the spectrum does indeed highly depend on 
the contributions of these extra partons. One observes a very good agreement between the 
NLL+NLO results and the Monte Carlo predictions after matching parton showering to 
matrix elements containing up to one extra parton, which is not surprising since the same 
matrix elements, relevant for hard emissions, are included in the resummation computations. 

The same type of effects can be noticed in Figures 8 and 9 for the benchmark scenarios 
18 and 31, the only difference being the overall scale of the distributions, which is about 
one order of magnitude smaller than in scenario 1 due to the larger gaugino masses. 

In Figures 10 and 11, we address the production of two other pairs of gauginos, i.e. of 
X]~X2 ^^^^ XiXi- As before, the Monte Carlo curves have been normalized according to 
the resummed results of Table 4, i.e., to 17.05 fb and 30.04 fb for the associated production 
and the chargino pair production channels, respectively. The distributions are similar to 
those of the golden channel but the absolute size of the latter (see Figure 7) exceeds 

the one of the xrX2 channel (see Figure 10) as expected for proton-proton collisions where 
the production of a positively charged final state is favored with respect to a negatively 
charged final state. 

4 Conclusion 

With a Higgs boson of about 125 GeV to be unveiled, but no hint yet for any colored 
supersymmetric particle, experimental searches are now more and more focusing on the 
pair production of the superpartners of the gauge and Higgs bosons. 

In this paper, we have updated the total cross sections related to the production of 
various gaugino pairs at the LHC, running presently at a center-of-mass energy of 8 TeV. 
We have presented results at the leading order and the next-to-leading order of perturbative 
QCD and matched, in addition, the NLO results to threshold resummation at the next- 
to-leading logarithmic accuracy in order to consistently account for the large logarithmic 
corrections arising from soft parton emission. Our cross sections have been given for several 
benchmark scenarios motivated by the recent LHC supersymmetric searches, together with 
the theoretical uncertainties stemming from scale variation and the choice of the parton 
density sets. We have shown that resummation at the NLL accuracy allows for a drastic 
reduction of the scale dependence of the cross section to less or about a percent, in contrast 
to the PDF errors which are now the dominant source of theoretical uncertainties and in 
general of the order of 2-3 percent. 

We have also analyzed the precision of the traditional experimental approach for su- 
persymmetric process simulation consisting in merging parton showering with multi-parton 
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Figure 8. Same as Figure 7 for the benchmark point 18. The gaugino masses are here 472 GeV. 



matrix elements. We have shown that the predictions obtained fohowing this approach are 
largely in agreement with the most precise theoretical computations matching resummation 
at the NLL level and fixed order calculations at the NLO accuracy. 

We finally emphasize that even if increasing experimental limits on the squark and 
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Figure 9. Same as Figure 7 for the benchmark point 31. The gaugino masses are here 479 GeV. 



gluino masses might exclude the investigated scenarios, the (lower) gaugino masses may 
still remain allowed so that our results could stay valid in the context of more general 
SUSY-breaking scenarios yielding gaugino masses of the same order. As already stated 
above, further numerical results are available from the authors upon request. 
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Figure 10. Same as Figure 7 for the production of a Xi X2 associated pair. 

A Total cross sections at y/s = 7 TeV for benchmark points 1, 18, and 31 

In this Appendix, we present in Tables 7, 8, and 9 the total production cross sections for 
the same combinations of neutralino and chargino pairs in the context of the benchmark 
points 1, 18 and 31 as in Section 2.2, but now for data taken earlier at the LHC at a 
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center-of-mass energy of = 7 TeV. We give again numerical results at the leading and 
next-to-leading order of perturbative QCD as well as after matching the NLO results to 
threshold resummation. 

Due to the slightly lower available energy, the cross sections are smaller by about 
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Process 



mi [GeV] 



ms [GeV] 



LO [fb] 



NLO [fb] 



NLO+NLL [fb] 



PP 
PP 
PP 
PP 
PP 
PP 
PP 
PP 
PP 
PP 
PP 
PP 
PP 
PP 
PP 



xixr 

x^x^ 

X2xr 

X3X4 
XsX^ 
X4X2' 
X^Xi 

X^X2 
X^X3 

x^xr 

X^Xs 

X^X4 

x^x^ 



161.7 
161.7 
303.8 
303.8 
303.8 
526.5 
526.5 
542.4 
303.5 
303.5 
303.5 
303.5 
542.2 
542.2 
542.2 



161.7 
303.5 
303.8 
526.5 
303.5 
542.4 
542.2 
542.2 
161.7 
303.8 
526.5 
303.5 
526.5 
542.4 
542.2 



61+^-^^" 
14+10.5% 

n cr9+12.2% 
"•"^^-10.2% 
n. 9c:+12.4% 
^•^^-10.4% 

n 94+12.5% 
"•^^-10.4% 
n 29+6.9% 

0.241-3? 

18.901-5 

n 01+12.4% 
^•^^-10.4% 

n 7f?+12.5% 
^•'"-10.4% 
Q t-4+12.1% 



0.79 
0.15 
0.73 
0.17 
12.33 
0.60 
0.30 
0.28 
0.35 
30.78 
0.28 
22.18 
0.91 
0.85 
0.62 



+3.7%+3.0% 
3.2%-2.1% 
+2.6%+3.1% 
2.6% -2.5% 
+3.7%+3.5% 
3.8% -2.2% 
+2.8%+3.6% 
3.2% -2.4% 
+1.7% +3.3% 
2.0% -2.7% 
+3.1% +4.4% 
3.6% -2.6% 
+3.1%+5.7% 
3.6% -4.2% 
+2.9%+5.8% 
3.5% -4.3% 
+2.6%+3.0% 
2.6% -2.2% 
+1.7% +3.2% 
1.9% -2.3% 
+2.8%+4.0% 
3.2% -2.6% 
+1.8%+2.9% 
2.1%-2.1% 
+3.1%+4.9% 
3.7% -3.0% 
+3.0% +4.9% 
3.6% -3.0% 
+2.9% +4.6% 
3.5% -2.7% 



+0.4% +2.9% 
-0.7% -2.1% 
n 14+0.0% +3.1% 
^•-^^-0.4% -2.7% 
+0.1%+3.7% 
0.1% -2.1% 
+0.2% +3.7% 
0.5% -2.4% 
+0.3% +3.3% 
0.8% -2.7% 
+0.6% +4.5% 
1.2% -2.8% 
+0.7% +5.8% 
1.3% -4.1% 
+0.8% +5.9% 
1.4% -4.2% 
+0.1% +3.0% 
0.4% -2.2% 
+0.0% +3.2% 
0.4% -2.3% 
+0.0% +4.1% 
0.3% -2.6% 
+0.1% +2.9% 
0.6% -2.1% 
+0.5% +4.8% 
-1.0% -3.0% 
+0.6% +4.9% 
1.1% -3.3% 
+0.6% +4.6% 
-1.2%-2.7% 



0.77 



0.71 

0.17 

12.18 

0.59 

0.30 

0.28 

0.35 

30.47 

0.27 

21.93 

0.91 

0.85 

0.61 



Table 7. Total cross sections related to the production of various gaugino pairs of masses mi 
and 1712, presented together with the associated scale and PDF uncertainties for the LHC running 
at a center-of-mass energy of y/s = 7 TeV in the context of the benchmark point 1 of the LPCC 
numbering scheme. The cross sections are given at the leading order and next-to-leading order of 
perturbative QCD and matched to threshold resummation. The PDF uncertainties are not shown 
for the LO results. Any cross section smaller than 0.1 fb is omitted. 



Process 


mi [GeV] 


m2 [GeV] 


LO [fb] 


NLO [fb] 


NLO+NLL [fb] 


PP X2X1 
PP xtX2 

PP xtxi 


471.9 
471.8 
471.8 


471.8 
471.9 
471.8 


9 n7+ii.o% 
^•^'-9.3% 


1 9i+2.0%+4.7% 
i.Zi_2 7%_3 5% 

q cq+2.0%+4.3% 
0.00_2 7%_2 7% 

9 qq+2.1%+3.9% 
Z.OO_2g%_2.4% 


1 -iQ+O.8%+4.8% 
J--J-y_i.4%_3.4% 

q 40+0.6% +4.4% 
"^•^^-1.1% -2.8% 
9 qn+0.6%+3.8% 
^•"5^-1.3%-2.5% 



Table 8. Same as Table 7 for the benchmark point 18 of the LPCC numbering scheme. 



30% with respect to those presented in Tables 4-6. In particular, the cross sections for 
heavy gaugino pairs do not exceed 0.1 fb at the benchmark points 18 and 31, so that they 
are no longer listed. On the other hand, the ATLAS and CMS experiments accumulated 
in 2011 an integrated luminosity of around 6 fb~^ each, so that for the light channels a 
significant number of gaugino events could have been produced at all three benchmark 
points considered here. 
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Process 


mi [GeV] 


ma [GeV] 


LO [fb] 


NLO [fb] 


NLO+NLL [fb] 


pv xlxi 

PP xtx2 

pp xtxi 


478.5 
478.5 
478.5 


478.5 
478.5 
478.5 


2 04+11-1% 
1 qn+11-1% 


1 -in+2.0%+4.8% 

i.iU_2,7%_3.5% 

q 90+2.0% +4.4% 
0.zo_2 7%_2,7% 

9 10+2.1% +3.9% 
Z.iO_2,8%_2.4% 


1 ns+O-9%+4-6% 
-^•^°-1.5%-3.6% 
9r,+0.4%+4.6% 
"^•^^-1.1% -2.8% 
9 -,1+0.5% +3.8% 
^•-^^-l.l%-2.6% 



Table 9. Same as Table 7 for the benchmark point 31 of the LPCC numbering scheme. 
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